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Abstract. Exact results are derived on the averaged dynamics of a class of random 
quantum-dynamical systems in continuous space. Each member of the class is characterized 
by a Hamiltonian which is the sum of two parts. While one part is deterministic, time- 
independent and quadratic, the Weyl-Wigner symbol of the other part is a homogeneous 
Gaussian random field which is delta correlated in time, but smoothly correlated in position 
and momentum. The averaged dynamics of the resulting white-noise system is shown to be 
a monotone mixing increasing quantum-dynamical semigroup. Its generator is computed 
explicitly. Typically, in the course of time the mean energy of such a system grows linearly 
to infinity. In the second part of the paper an extended model is studied, which, in 
addition, accounts for dissipation by coupling the white-noise system linearly to a quantum- 
mechanical harmonic heat bath. It is demonstrated that, under suitable assumptions on 
the spectral density of the heat bath, the mean energy then saturates for long times. 

Keywords. Random quantum-dynamical systems, white noise, dissipation. 

1 Introduction and summary of results 

Ever since Anderson (1958) one-particle quantum Hamiltonians with time- 
independent random potentials have been intensively studied as models for 
disordered electronic systems (Bonch-Bruevich et al., 1984; Shklovskii and 
Efros, 1984; Lifshits et al., 1988; Carmona and Lacroix, 1990; Pastur and 
Figotin, 1992). In contrast, models with time-dependent random potentials 
have attracted less attention. They have been used, for example, to explain 
the diffusive behaviour of Frenkel excitons in molecular crystals (Haken and 
Rcineker, 1972). Subsequently Ovchinnikov and Erikhman (1974), Madhukar 
and Post (1977) and Girvin and Mahan (1979) studied lattice models with 
time-dependent diagonal and off-diagonal random parts. Here the random- 
ness enters via time-dependent fluctuating coupling coefficients and serves to 
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mimic the influence of collective excitations of the crystal. At high tempera- 
tures the coupling coefficients should be uncorrelated in time. More precisely, 
it is justified (Ovchinnikov and Erikhman, 1974) in this case to view them 
as realizations of a Gaussian white noise. One is then able to derive exact 
analytical results on the averaged dynamics of the model. 

The simplifying feature of a vanishing correlation time was also exploited 
by Jayannavar and Kumar (1982). They were concerned with the averaged 
dynamics of a point particle moving in continuous space under the influence of 
a time-dependent Gaussian random potential which is uncorrelated in time, 
but smoothly correlated in position. Their main findings are a cubic long- 
time growth for the spatial spreading of states under the averaged dynamics 
and a linear increase in time of the mean energy. This behaviour is strikingly 
different from the diffusive spreading of states known from the aforementioned 
lattice models, which goes along with a saturating mean energy and may be 
interpreted in terms of Umklapp-pr ocesses. An extension of Jayannavar and 
Kumar (1982) to the case of an additionally present constant electric field is 
given in Jauho (1987). 

A detailed study of models with uncorrelated time-dependent random- 
ness is not only of interest for understanding specific problems in condensed- 
matter physics. Indeed, they are of relevance to many different branches of 
physics. Accordingly, the appearance of an additive white noise is not nec- 
essarily due to thermal fluctuations. Such models occur, for instance, in the 
debate on fundamental issues of quantum mechanics. In this context their 
averaged dynamics allows for a quantum-mechanical description of classi- 
cal objects while still forbidding superpositions of macroscopically distinct 
state vectors, see for example Joos and Zeh (1985), Ghirardi et al. (1986), 
Joos (1987), Ghirardi et al. (1987), Bell (1987, p. 201), Ghirardi and Rimini 
(1990), Zurek (1991), Pearle and Squires (1994). In the realm of cosmology 
they are sought (Hawking, 1982; Banks et al., 1984; Hawking, 1984) to solve 
the information problem (Hawking, 1976) which arises by the evaporation of 
black holes. Furthermore they are met in classical kinetic theories (Spohn, 
1980) and in quasi-classical Markovian transport theories (Markowich et al., 
1990; Ferry and Grubin, 1995; Haug and Jauho, 1996). 

Hence it is desirable to provide a systematic discussion of the stochastic 
dynamics of a rather general model with uncorrelated time-dependent ran- 
domness. This was initiated in our earlier paper (Fischer et al., 1994), where, 
as a first step, the averaged dynamics of such a model was studied. Since 
proofs of the results were only sketched there, it is one goal of the present 
paper to provide detailed proofs. The class of quantum-dynamical random 
systems which was introduced by Fischer et al. (1994) — see Model 2.1 in the 
next section — is characterized by a model Hamiltonian in continuous space 
which is the sum of a deterministic part and of a random part. Apart from 
allowing the deterministic part to be quadratic in both momentum and posi- 
tion, which may serve to describe the influence of a constant magnetic field, 
say, it is the other main feature of Fischer et al. (1994) that the random part 
may have fluctuations in both position and momentum. The fluctuations 
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in momentum may be viewed as a continuum analogue of the off-diagonal 
random parts in the lattice models mentioned above or as a caricature for 
inelastic collisions of the particle. From the point of view of canonical me- 
chanics they are natural anyway. 

As to the mathematical description of random operators with fluctuations 
in position and momentum, one could think of applying a factor-ordering 
prescription to a random field on classical phase space. Therefore, and in 
order to clearly isolate quantum effects from classical ones even before taking 
expectation values, it is only natural not to use the conventional Hilbcrt-space 
formulation of quantum mechanics, but to represent the quantum system 
from the outset in phase space. In so doing, we choose the linear phase- 
space representation of Weyl, Wigner and Moyal, which amounts to a totally 
symmetric operator ordering when going back to Hilbcrt space. 

As was already mentioned it is one goal of the present paper to provide 
detailed proofs of the results given in Fischer et al. (1994). This is done 
in Section 2. In Statement 2.2 we establish the semigroup property of the 
averaged dynamics of Model 2.1. Its irreversible character is described in 
Statement 2.4. Moreover, the averaged dynamics is positivity preserving and 
hence even a Markov semigroup, see Remark 2.6. Statement 2.7 is concerned 
with some geometrical properties of averaged time-evolved states and State- 
ment 2.11 shows that in typical cases the mean energy of Model 2.1 grows 
linearly to infinity in the course of time. This unlimited increase of energy oc- 
curs because there is no mechanism built into Model 2.1 which causes energy 
to dissipate. It is therefore the second goal of the present paper to study also 
an extended version of Model 2.1 in which dissipation is incorporated, too. 
This is done in Section 3, where the white-noise system is linearly coupled 
to a quantum-mechanical heat bath of independent harmonic oscillators, see 
Model 3.1. The averaged dynamics of the total system, that is, white-noise 
system plus heat bath, is analyzed in Subsection 3.1 with the methods devel- 
oped in Section 2. In Subsection 3.2 we eliminate the degrees of freedom of 
the heat bath assuming a thermal equilibrium state of the heat bath. State- 
ment 3.12 at the end of Subsection 3.4, which is the main result of the second 
part of the paper, shows that the thus obtained reduced averaged dynamics 
of Model 3.1 acquires the desired dissipative features in the macroscopic limit 
of the heat bath provided its spectral density behaves suitably. The interplay 
of noise and dissipation then prevents the mean energy of the white-noise sys- 
tem from diverging in the long-time limit. The proof of this result requires 
some technical preparations which are dealt with in Subsections 3.3 and 3.4. 
Section 4 compiles some references concerning two related challenging prob- 
lems which we do not address otherwise in this paper: fluctuations around 
the averaged dynamics and effects of a non-vanishing correlation time, that 
is, coloured-noise perturbations. Finally, for the convenience of the reader, 
the basic properties of the Weyl-Wigner-Moyal representation of quantum 
mechanics are presented in the Appendix. The reader who is interested in 
more details is referred to the (hopefully) up-to-date list of references. 
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2 The averaged dynamics by white-noise Hamiltonians 

In what follows, we are interested in the non-relativistic quantum dynamics 
of a single spinless point particle whose configuration space we choose to be, 
for simplicity, the Euclidean line R . Following Weyl, Wigner and Moyal we 
will represent this quantum system in the two-dimensional symplectic phase 
space IxR. The generalization to several Cartesian degrees of freedom is 
straightforward and will be dealt with on demand in Section 3. 

2.1 White-noise Hamiltonians 

In this subsection we define the basic model of this paper. 

Model 2.1. Being the Weyl-Wigner symbol of a Hamiltonian operator, the 
Hamiltonian function 

H + N t , (2.1) 

also referred to as white-noise Hamiltonian in the sequel, defines a random 
quantum- dynamical system on phase space IxM. The deterministic part 
H of (2.1) is time independent and at most quadratic in momentum p and 
position q. The time- dependent random part N t is a homogeneous Gaussian 
random field in p, q and time t with zero mean, 



N t (p,q)=0, (2.2) 
and covariance function 



N t (p,q)N t ,( p i,qi)=6(t-t')C(p-p / ,q-q'). (2.3) 

Here we denote the stochastic average over randomness by an overbar, d 
stands for the Dirac delta function and the phase-space part C of the covari- 
ance function is assumed to be arbitrarily often differentiable. 

The requirement for H to be at most quadratic in position and momentum is 
necessary for deriving explicit results on the averaged dynamics of Model 2.1. 
On the other hand the time independence of H and the homogeneity of Nt are 
only assumed for simplicity. If one were interested in treating multiplicative 
noise (Habib, 1994), this could be managed, for example, by giving up the 
homogeneity assumption. 

We would like to stress that in order for C to be a covariance function it 
needs to be the Fourier transform of a positive measure on 1 x R . This point 
is disregarded by Jayannavar (1993). Therefore it comes without surprise 
that some of his results are physically insignificant — especially those which 
are in contradiction to ours. 

The dynamics generated by the white-noise Hamiltonian (2.1) is formally 
described by the stochastic quantum Liouville equation (Kubo , 1963) 



d t w t = 



-[H,w t ]-[N u w t ] 



(2.4) 
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with initial condition wq = w. Here w is the Wigner density representing 
the initial state and [ • , • ] denotes the Moyal bracket (A7) of two phase- 
space functions. Note that the first Moyal bracket [H,wt] in (2.4) is in fact 
a Poisson bracket due to the quadratic nature of H, cf. Remark A. 3. Being 
a stochastic partial differential equation, (2.4) has to be supplemented by a 
prescription of how to interpret the second Moyal bracket on its right-hand 
side, which involves the distribution- valued white noise N t . We choose a 
prescription which is quite natural on physical grounds: The time change 
of a realistic physical system is governed by driving forces with a non-zero 
correlation time. Theoretical models based on stochastic processes which are 
uncorrelated in time may only be used successfully if the correlation time of 
the actual driving forces is much smaller than any other time scale inherent 
to the system. Therefore we interpret (2.4) in the Stratonovich sense which — 
as is supported by the Wong-Zakai-like theorems in Horsthemke and Lefever 
(1984, p. 101), Karatzas and Shreve (1988, Chap. 5.2.D) and Brzezniak and 
Flandoli (1995) — corresponds to considering the white-noise limit t\0 of 

(t) 

the solution w t of the regularized initial-value problem 



dM T) =-[H,w^]-[Nt\w[ T) ] 



(2.5) 



In (2.5) we have replaced the Gaussian white noise Nt by the Gaussian 

coloured noise N^ which vanishes in mean, Nj: T \p,q) = 0, and is corre- 
lated according to 



iV t (T) (p, q) N[r\p', q') = g {T \t - t') C(p -p',q- q') , (2.6) 

where is a correlation function that approaches the Dirac delta function 
in the limit of vanishing correlation time r\0. For times t > the solution 
of (2.5) can be conveniently expressed in terms of a time-ordered exponential 

w { t T) =w + yV-l) n / ds n / da n _i---/ dsi 
n=1 Jo Jo Jo 

x[H + ATM [H + N^,. ..,[H + N£\w] . . .]] 
= :cxpj-^ ds [H + N ( /\ .]|«7. (2.7) 

We remark that the corresponding solution for times t < requires the time 
ordering to be reversed. 

2.2 The averaged dynamics 

The averaged dynamics generated by a white-noise Hamiltonian is the map- 
ping % defined by 

w^T t w:~wl, (2.8) 
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where wi(p, q) := w t (p, q) denotes the stochastic mean of the solution of (2.4) 
in the Stratonovich interpretation. According to the preceding subsection, T t 
can be expressed as the limit 



% = lim exp|-y ds [H + 7Vj r) , ,]J , t>0. (2.9) 

The salient feature of the white-noise limit is that it allows to describe the 
averaged dynamics explicitly. 

Statement 2.2. The averaged dynamics of H + N t is a semigroup of linear 
operators 

% =cxp{-i(£+A0}. (2.10) 
Its generator consists of the Liouvillian 

C:=[H, .} = (d p H)d q - (d q H)d p (2.11) 

associated with the deterministic part H of the Hamiltonian and of the per- 
turbation 

N := h- 2 {C(0,0) - C(-iHd q ,ihd p )} (2.12) 
which depends only on the stochastic properties of the random part N t . 

Proof. To perform the stochastic average in (2.9) it is advantageous (Kubo, 
1962, §7 (5)) to isolate the unperturbed time evolution generated by C in 
passing to the Dirac-Dyson representation of the time-evolution operator 

expj-^ ds (£ + [iVj T) , .])| =e~* £ expj-^ ds e s£ [iVj r) , .]e- s£ |. 

(2.13) 

Next we introduce the family of Gaussian-distributed random complex mea- 
sures Nj; T \dpdq) on phase space arising in the spectral representation (Doob, 
1964, Chap. 11, § 4) 



(r) 



[p,q)=: f N^ ) (dp'dq')c i ^ q '- q P'^ h (2.14) 

JKxI 



of N^ T \ They have zero mean, Nj: T \dpdq) — 0, and the covariance measure 



Ni T> (dpdq) N^'idp'dq') = gW(t - t') C{dpdq) S_ P} _ q (dp'dq') . (2.15) 

Here S- p .- q (dp'dq') = 5{p' + p)6(q' + q)dp'dq' denotes the Dirac measure 
with support (— p, —q) elxl. The covariance function C appears as the 
symplectic Fourier transform 

C(p,q)= [ C(dp'dq') jW-qp')/* (2.16) 
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of the measure C on R x R . Inserting (2.14) into the definition (A7) of the 
Moyal bracket one obtains 

[N ( /\ .) = i / N^(dp'dq') F(p',q') . (2.17) 

JRxK 

Here we have introduced the operator-valued function 

lxR3 (p', <?') -> ^(p', 9') := h- 1 e i W-oi>'» h (V p , /2 . q//2 - V_ p , /2 ,_ q>/2 ) . 

(2.18) 

In (2.18) p and q play the role of the canonical coordinate functions on phase 
space and V p ' l9 / := e p dp+q dq is the translation operator by (p 1 , q'). In other 
words, T(jl ', q') acts on a phase-space function / according to 

{T(p\ q')f) (p, q) = h- 1 jW-rfV* {f(p + p'/2, q + q> /2) 

-f(p-p'/2,q-q'/2)}. (2.19) 

With these preparations and writing J r s (p',q l ) := c sC J r (p',q')e~ sC for the 
Dirac picture of !F{p' , q'), the averaged dynamics can be cast into the form 



Tt=e-* £ lim exp(-i / ds / N^\dp'dq') F s (p', q')\ . (2.20) 

rX0 I JO JlxK J 

To perform the stochastic average in (2.20) we appeal to the characteristic 
functional of Gaussian random variables and conclude that 

T t = c- tc Urn cxp j-i J ds ds' g {r \s - s') 

x f C(dp'dq') F s (-p',-q') F s ,(p',q')\ 

JrxR J 

exp(-i fds [ C(dp'dq')T s (-p',-q')T s (p',q')\. (2.21) 
I 1 JO ilxl J 



Furthermore one has 



/ C(dp'dq')f(-p',-q')F(p',q') 

JRxR 

= h- 2 f C(dp'dq') (2 - V pW - V-j,,-*) 
ilxl 



= 2h- 2 (C(0,0) - C{-\hd qi ihd p )} 

= 2Af, (2.22) 



because C is an even function. Hence 



T t = e-* £ cxpj-^ ds c s£ AAc- s£ | = 0^+^ , (2.23) 



and the proof is complete. 
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Remark 2.3. It follows from Statement 2.2 by a differentiation with respect to 
time that the averaged state wi obeys the linear integro-differential equation 

dm(p, q) = [m, H](p, q) + H- 2 I C(dp'dq') {wi(p + p' , q + q') 

JRxR 

-Wi(p,q)} (2.24) 

with initial condition Wo — w. Eq. (2.24) is a generalization of the main result 
of Jayannavar and Kumar (1982). In the special case of a free deterministic 
part, H = p 2 /2m, and a momentum-independent noise, that is, when C does 
not depend on p, it reduces to an equation which is equivalent to Eq. (8) in 
Jayannavar and Kumar (1982). The latter is derived there by directly aver- 
aging the von Neumann equation with the help of a partial integration with 
respect to the Gaussian probability measure (Furutsu, 1963; Donsker, 1964; 
Novikov, 1964), see also Rytov et al. (1989, p. 32) or Glimm and Jaffc (1987, 
Thm. 6.3.1). The present approach is more direct and appears to provide a 
better starting point towards a generalization to the case of non-vanishing 
correlation times, that is, coloured noise. Note also that the subsequent 
treatment of Eq. (8) in Jayannavar and Kumar (1982) is restricted to a pure 
initial state which, in Weyl-Wigner language, is represented by a Gaussian 
Wigncr density w. 

Being the stochastic average of a Hamiltonian dynamics the semigroup % 
is completely positive and thus provides an example of a quantum-dynamical 
semigroup in the sense of Bausch (1966), Gorini ct al. (1976), Lindblad 
(1976), Spohn (1980), Alicki and Lendi (1987) or Streater (1995). Quantum- 
dynamical equations similar to (2.24) are discussed in very different branches 
of physics. One reason for this is revealed by 

Statement 2.4. The averaged dynamics is monotone mixing increasing in 
the sense that 

d t (T t w,T t w) <0 (2.25) 

for allt > and all Wigner densities w. Here we use the notation ( . , . ) for 
the standard scalar product (A4) of phase-space functions. In (2.25) equality 
holds if and only if C is a constant. 

Proof. Taking the scalar product of (2.24) with Wi and observing 
(Wi, [wi, H]) = 0, which follows from a partial integration, we get 

(Wi, d t Wi) = h- 2 [ C(dp'dq') {(Wi, V p >, q >Wi) - (wi,Wi)} ■ (2.26) 

Due to the positivity of the measure C, inequality (2.25) now follows from 
an application of the Cauchy-Schwarz inequality 

(wi,V p >, q >wi) < (wi,Wi) 1/2 (Vp>, q >wi,V p >, q >wi) 1/2 = (Wi,wi) . (2.27) 
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To obtain equality in (2.25) it is necessary and sufficient that there is equality 
in (2.27) for (7-almost all (p' 7 q r ). Hence, it is necessary and sufficient that 
the functions wi and V p ' iq >wi are linearly dependent for C-almost all (p', q'). 
In this case there are constants \ p ' , q > with V p iywi = \ p >. q 'Wt- Normalization 
{l,V p ' t q'Wt) = (l,iOt) requires \ p >, q ' = 1. For (p',q') ^ (0,0) this implies 
periodicity of the Wigner density wi which is in contradiction to its square- 
integrability (wi,wi) < < oo, cf. Statement A. 5. To summarize, 

the Wigner densities wi and V p ' iq 'W~i are linearly dependent if and only if 
(p'> ?') — (0, 0) which in turn allows for equality in (2.25) if and only if C has 
support {(0, 0)}. ■ 

Statement 2.4 can be reformulated as dtS2{%w) > 0, where 52 (w) := 
— ln(27r/i (w, to)) is the Renyi entropy of order two of the quantum state 
(Jumarie, 1990). Similar to the von Neumann entropy, S2 reflects the degree 
of mixing of the quantum state (Thirring, 1983, (2.2,3)) and, by Jensen's 
inequality, it is a non-negative lower bound to the former. The average over 
randomness has turned the fully reversible quantum Liouville equation (2.4) 
into the "master equation" (2.24) with coherence-destructing irreversible be- 
haviour. 

So far we have not made use of the fact that the deterministic part H of 
the Hamiltonian is at most quadratic in p and q. The next result however 
relies crucially on it. 

Statement 2.5. The averaged dynamics admits the representation 
rt 



% = e~ tc 



cxp j-fr 2 J ds ((7(0, 0) - C{-\hK s , ihX s )) I . (2.28) 



The coefficients of the first-order differential operators 

jc s := (d q c- sC P )d p + {d qC - sC q)d q , 

X, : (d p c- sC p)d p + (O p e- sC q)d q 

depend only on time and not on position and momentum. They are deter- 
mined by the classical phase-space trajectories generated by H. 

Note that there is no time ordering in (2.28) because the four operators /C s , 
K s i, X s , X s i commute pairwise due to the quadratic nature of H. Statement 
2.5 is important because it allows to derive explicit results on the averaged 
dynamics %. 

Proof of Statement 2.5. Thanks to (2.23) it only remains to show that 

e sC Afe- sC = h- 2 (C(0, 0) - C(-itUC a ,ihX„)) . (2.30) 

But when applying both sides of (2.30) to a phase-space function, this identity 
follows from the chain rule and from observing that the phase-space trajec- 
tories {(e* p)(p' , q'), (e q)(p', q')} t>0 are linear functions in the initial data 
(p' } q') due to the quadratic nature of H. m 
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To gain some more insight into the behaviour of % define the phase-space 
function 

P t (p,q) : = e -tc(o,o)/Hi exp j ft - 2 ^ ds (Cojt_ s )( p , q) J , (2 . 3 1) 
where 

Jf-=( d fl C c P d /tl P ) (2-32) 
1 o p e tL -q o q e tL -q 1 



is 



is the Jacobian of the flow e tc on phase space and -4 := ( i V ) ^ ( -i o ) 
its adjoint with respect to the canonical symplectic form ((p,q),(p' ,q')) i— > 
P<z' — QP 1 - The superscript T denotes the transposition of matrices and the 
circle o denotes the composition of mappings. Because of 



— ihK, s \ ,t f—\hd q 



\hX s I s \ ihd p 



(2.33) 



one can write % — e~ tc P t (—ihd q ,ihd p ). Now observe that P± is a charac- 
teristic function, that is, the Fourier transform of a measure P t on phase 
space 

P t (p,q)=: [ P t (dpW)c i(p "'- w ' )/R . (2.34) 

This is true because h~ 2 J *ds CoJ^ s is a characteristic function and because 
sums and products of characteristic functions are again characteristic func- 
tions. In fact, P t is even a probability measure due to P t (0,0) = 1. Thus we 
conclude 

% = e- tc [ Pt(dp'dq r ) V pl , ql . (2.35) 

Remark 2.6. The result (2.35) shows that the averaged dynamics is also pos- 
itivity preserving and thus even a Markov semigroup. In particular, it maps 
classical states, that is, probability densities onlxl, to classical states. 
Eq. (2.24) appears as the associated kinetic equation, which — up to the drift 
term arising from H — may be viewed upon as a linear Boltzmann equation 
with a homogeneous stochastic kernel (Lasota and Mackey, 1994; Reed and 
Simon, 1979, Chap. XI.12). 

Since the covariance function C of the random field is an even function, all 
odd moments of the probability measure Pt vanish. The covariance matrix 
of P t 



s 
s 

Jo 



P,(dpd9)(£ 5)=-» a ((fl a P«)(0,0)) t 
dsJ_ s ((-9 2 C)(0,0)) f JT S (2.36) 
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is monotone increasing in time, 

0<C t /<Q, 0<t'<t. (2.37) 

In (2.36) we have used the symbol d 2 for the Hessian of second-order partial 
derivatives. For later purpose it is convenient to introduce the non-negative 
matrix of diffusion coefficients 



whose entries 



:= {(-id P Y{id q yc) (o, o)/( M ! v \) (2.39) 

are given by the curvature of the covariance function at the origin. Thus one 
can also write 



i 



t 

daJ-.DJl,. (2.40) 



It is clear from the Riemann-Lebesgue lemma that the probability mea- 
sure Pt cannot be absolutely continuous with respect to Lebesgue measure 
dpdq. The following decomposition helps to decide whether P t possesses an 
absolutely continuous component. 

Statement 2.7. The averaged dynamics of an initial state w 
(T t w)(p,q) = e-* c (°'°)/ ft2 w(e- tc p,e- tc q) 

+ [ Qt(dp'dq') w(p' + e- tc p,q / + C~ tC q) (2.41) 

decomposes into an exponentially decaying contribution of the unperturbed 
dynamics and into a contribution where the initial state is first smeared out 
by the finite measure Q t and then subjected to the unperturbed dynamics. The 
measure Q t is defined by its characteristic function 

Q t := P t - e -* c (°'°)/ ft2 (2.42) 

and is normalized to Q t (0,0) = 1 — c - tc '(v<°)/ h2 . All higher moments of Q t 
coincide with those of P t . A sufficient condition for Q t to have an integrable 
and continuous density with respect to Lebesgue measure is the integrability 
of C. The density is then bounded by 

°<^^<7^W dp'dq' \C(p',q')\. (2.43) 
dpdq (2irh 2 ) 2 7 RxR 

Necessary conditions for the existence of a density are the strict positivity of 
the covariance matrix, Q > 7 and the vanishing of the function (p, q) 
J *ds (CoJ^ s )(p, q) at infinity. 
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t = t2>tl 










J L° 


A / 










t = ti > 



Fig. 1. Illustration of the averaged dynamics of a Wigncr density in phase space. The 
deterministic part of the Hamiltonian has been chosen to be that of a harmonic oscillator. 
The solid (dashed) contribution symbolizes the first (second) term in the decomposition 
(2.41). 

Figure 1 illustrates the situation described in Statement 2.7 for the case where 
H is the Hamiltonian of a harmonic oscillator and w a non-negative initial 
Wigner density at time t = 0. The dotted ellipse represents one phase-space 
trajectory of H. The averaged state %w at times t\ > and t 2 > t\ is ob- 
tained by summing up the solid and the dashed contribution which symbolize 
the first and second term in the decomposition (2.41), respectively. Accord- 
ingly, the weight of the solid contribution to %w decreases exponentially 
in time, whereas that of the dashed one increases to one. In addition, the 
width of the dashed contribution increases in time. The sum of the volumes 
enclosed underneath is always one. 

Proof of Statement 2.7. The arguments which established that P t is a char- 
acteristic function apply to Qt as well. Hence, Qt is a finite measure on 
Ixl. According to Bauer (1996, § 23, Problem 10) the integrability of Q t 
is sufficient to guarantee the existence and continuity of the Lebesgue den- 
sity Q t (dpdq)/dpdq. But the inequality \e z — 1| < \z\e y , which is valid for 
all \z\ < y, shows that this condition is fulfilled if C is integrable. Here, we 
have also used the fact that phase-space volume is conserved by Hamiltonian 
flows. The upper bound in (2.43) then follows from the Fourier inversion for- 
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mula. The second of the two necessary conditions is just the content of the 
Ricmann-Lebesgue lemma. As to the first one we remark that if it did not 
hold, there would exist a pair [p' , q 1 ) such that J RxR Q t (dpdq) (p'q—q'p) 2 = 0. 
But this implies that Q t is supported on a set of Lebesgue measure zero. ■ 

2.3 The classical limit 

One main advantage of a phase-space representation of quantum mechanics 
is that the classical limit of a quantum system may be studied with little 
efforts. We denote the averaged dynamics of the classical white-noise system 

by 

r t , c i := limT t (2.44) 
and conclude 

Statement 2.8. Suppose that the deterministic part H of the Hamiltonian 

and the covariance function C do not depend on h~. 

Then the averaged dynamics of the classical white-noise system 

i) is a semigroup of linear operators 

T tiC i =exp{-i(£+A/;i)} = e- t£ expji (d p ,d q ) Q (J) j (2.45) 
with generator C + 7V c i , where 

Ki~-(d P ,d q )DQ , y (2.46) 

The covariance matrix Q and the matrix D of diffusion coefficients were 
introduced in (2.36) and (2.38), respectively. 

ii) maps an initial probability density p on phase space to the probability 
density ~pi := Tt,c\ P which obeys the Fokker-Planck equation 

dm= ((d q H)d p - (d p H)d q )Tt + (d P ,d q ) D Q^jjH (2.47) 

with initial condition po = P- The solution of (2.47) is obtained by 
smearing out p with the Gaussian probability measure P tyC \ on phase 
space, which is defined by its characteristic function 

exp{-i (x, k) Q ( X \ j =: f P t>cl (dp'dg') e^*^') , (2.48) 

and subjecting the result to the time evolution of H 

Tt(p, q)= I PtA&P'M') Pip' + c- tc P , q' + c- tc q) . (2.49) 

JRxR 
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The measure P tyC \ has zero mean and C t as its covariance matrix. It 
possesses the Lebesgue density 

5^21 . (dct^C,))-" 2 expj-lMC- (>)} (2.50, 

if and only if Q > 0. 

Hi) is irreversible, if and only if the covariance function C is non-constant. 
In this case the Boltzmann-Gibbs entropy 

S(Tt) := -(TtM{2nKfH)) (2.51) 
of~pi (if existent) is a strictly increasing function in time. 

Proof. The first part of Statement 2.8 follows from Taylor expanding Af in 
(2.23) with respect to h 

M = -D , 2 d 2 p - D hl d p d q - D 2 , d 2 q + 0{h 2 d A ) , (2.52) 

with Dftv being defined in (2.39). In the limit fi\0 the expansion (2.52) 
reduces to the semielliptic differential operator M c \. To obtain the right 
equality in (2.45) it is useful to note in addition the operator identity 

dp ^e~ sC = C - sC lI s ( d A, (2.53) 



dj ~ s \d q/ 

which is a consequence of the chain rule and may be proved by applying both 
sides of (2.53) to a phase-space function. The second part of Statement 2.8 is 
an immediate consequence of the first part. To prove the third part, suppose 
that the Boltzmann-Gibbs entropy S{T t . c \ p) exists. Since T t ^ c \ is a semigroup 
it suffices to show S(%^ c \p) > S(p) if t > and C is non-constant. To do so 
we conclude from (2.49) 

S{T t>d p) = - ( P t . cl (dp'dq') (p, ln(2irhV- p ,,- q ,c tc T t , cl p) ) 

>-[ Pt.Mdp'dq') (pM^hp)) 

JRxR 

= S(p) , (2.54) 

where V p ' i(? < = e p dp+q dq denotes the translation operator by (p',q'). The 
strict inequality in (2.54) follows from the elementary inequality 

u — u\nu < v — u\nv , u,v>0,u^v. (2.55) 

To see this we note that for t > and a non-constant C, the support of 
P t , c \ is at least a one-dimensional subspace ofRxR. The integrability of p 
therefore forbids the equality 

p(p, q)= [ P t . cl (dp"dq") p(p" -p'+p, q" - q' + q) (2.56) 

JrxR 
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for pdpdq- almost all (jp,q) £ Ixl and P t>c i(dp'dg')-almost all (p',q') £ 
R x E. ■ 

Remark 2.9. The classical phase-space trajectories (pt,qt) generated by the 
Hamiltonian H + N t realize a diffusion process on phase space. Its underlying 
stochastic differential equation coincides with Hamilton's equations of motion 

d tPt = -(d q (H + N t ))(p t ,q t ), 

d t qt = {d p (H + N t ))(p t ,q t ), 

and (2.47) is the associated Fokker-Planck equation. 

2.4 Explicit results for averaged expectation values 

Now we return to the quantum-mechanical situation. In order to compute 
the averaged expectation value of an observable a, it is often advantageous 
to switch to the Heisenberg picture according to 

(w t , a) = {Wi, a) = (T t w, a) = (w, T t f a) . (2.58) 

The thus defined adjoint semigroup is obtained from T t by simply revers- 
ing the sign of C. In analogy to Statement 2.5 the adjoint semigroup has the 
representation 

7^ = e t£ exp|- I ds e~ sC N e sC \ (2.59a) 
= e tc 



o 



cxpj-^ 2 J ds (C(0, 0) - C(-ifi/C_ s , ihX- s )) j . (2.59b) 

Statement 2.7 finds a similar translation. The expression (2.59b) serves as 
our starting point for obtaining some explicit results on the averaged dynam- 
ics. For example, the averaged Heisenberg picture of any observable which is 
a polynomial in p and q may be calculated from it. To this end we remark 
that there are only finitely many terms in the Taylor expansions of the right 
exponential in (2.59b) and of the covariance function which contribute to the 
result. Examples for the cases where the deterministic part of the Hamilto- 
nian is that of a free particle or that of a particle subjected to a constant 
magnetic field in two space dimensions have been reported in Fischer et al. 
(1994). Here we only recall the free-particle case. 
Example 2.10. Let H = p 2 /(2m). Then 

T t t p = p, (2.60a) 
T t ] q = q + tp/m , (2.60b) 
T/p 2 = P 2 + 2tD aa , (2.60c) 
T t ] pq = p(q + tp/m) + tD 1A + t 2 D , 2 /m , (2.60d) 
T t V = (q + tp/m) 2 + 2tD 2fi + t 2 D 1A /m + § t 3 D , 2 /m 2 , (2.60e) 
7;Y = (q + tp/m) 4 + I2{(q + tp/m) 2 d 2 (t) + 2h 2 d 4 (t) + (d 2 (t)) 2 } . (2.60f) 
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The coefficients D^ v were denned in (2.39) and we have introduced the poly- 

nomials d a (t) := > - — 11 v ' v — t u+1 of maximal degree a + 1. As a con- 
M/ ^{u + l)m v 

sequence, the long-time spreading of states under the averaged dynamics is 
diffusive in momentum and superballistic in position 

lim \ ((w- t ,p 2 ) - (w-t,p) 2 ) = 2A).2 , (2.61a) 
lim 1 ({W U q 2 ) - {W t , q) 2 ) = ^ . (2.61b) 

Example 2.10 illustrates important features valid for general quadratic H: 
Observables which are linear in p and q are not affected by the white noise N t . 
Moreover, noise-induced terms in averaged expectation values of quadratic 
observables are independent of the initial state. Assuming that the deter- 
ministic part H of the Hamiltonian and the covariance function C are h- 
independent, noise-induced effects affected by quantum fluctuations occur 
only for observables of at least fourth order in p or q. However, very special 
situations are needed for quantum effects to show up in the leading term for 
long times. For example, taking the observables p n or q n , 4 < n integer, one 
must require the phase-space trajectories of H to grow exponentially in time. 
These results generalize some earlier ones in Jayannavar and Kumar (1982), 
Hcinrichs (1983b), Ghirardi et al. (1986), Heinrichs (1992) and Jayannavar 
(1993). 

Due to the time dependence of the random part N t of the Hamiltonian, one 
does not expect energy conservation for white-noise systems. The observable 

H can be interpreted as the mean energy of the system at time t. For the 
special case considered in Example 2.10, the mean energy grows linearly in 
time, see Eq. (2.60c). This is also true under more general circumstances. 

Statement 2.11. For an at most quadratic deterministic part H, the mean 
energy changes at a constant rate: 

Tt H = H + t (D .2(d 2 H) + D hl (d p d q H) + D 2 . (d 2 H)) . (2.62) 

The rate does not allow for a dependence on the initial state. 

Proof. The claim follows from (2.59a), from e sC H = H, from 

-AfH = D oa {d 2 p H) + D hl (d p d q H) + D 2fi (d 2 q H) (2.63) 

and from N 2 H = 0. Note also that the right-hand side in (2.63) does not 
depend on p and q. ■ 

Unless H is unbounded from below, such as H = (p 2 — Dq 2 )/(2m) for D > 
D 0,2/^2,0) the mean energy increases to infinity in the course of time. To 
compensate for this effect and, if possible, to allow for an eventual approach 
to a stationary equilibrium state, dissipation has to be incorporated. This 
will be done in the next section by coupling the white-noise system to a heat 
bath. 
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3 White noise and dissipation 

This section is devoted to the interplay of white noise and dissipation. Such 
an investigation is physically sensible, because the additive white noise in 
Model 2.1 is in general not due to thermal fluctuations. We will show that 
the incorporation of dissipation causes the mean energy of the white-noise 
system to saturate for long times. Topics of this kind were already addressed 
to in Heinrichs (1983a), where dissipation is generated through an ad-hoc 
time dependence of the deterministic part of the Hamiltonian. In contrast, 
we introduce dissipation by letting the white-noise system interact with a 
microscopic model of a heat bath — a more physical strategy which is com- 
monly accepted today (Magalinskii, 1959; Feynman and Vernon, 1963; Ford 
ct al., 1965; Ullersma, 1966; Arai, 1981; Caldeira and Leggett, 1983; Haake 
and Rcibold, 1985; Ford et al., 1988a; Grabert et al., 1988; Sogo and Fuji- 
moto, 1990; Weiss, 1993; Efimov and von Waldenfels, 1994; Jaksic and Pillet, 
1997a). 

The programme of this section will be accomplished in four steps: First, 
the white-noise system is bilinearly coupled to a heat bath of independent 
harmonic oscillators. The averaged dynamics of the total system, that is, of 
the white-noise system plus heat bath, is then studied with the methods of 
the preceding section. Second, the bath variables are eliminated assuming 
that the bath is in thermal equilibrium. Third, the thus obtained reduced 
averaged dynamics is considered in the macroscopic limit of the heat bath. 
Finally, the long-time limit is analyzed. Statement 3.12 summarizes the main 
result. 

3.1 The averaged dynamics of the total system 

For notational simplicity, we will no longer allow for a momentum dependence 
of the white noise and we will restrict the deterministic part of the Hamil- 
tonian to that of a harmonic oscillator. Note also that in this section we 
are using the Weyl-Wigner-Moyal representation for a quantum-mechanical 
system with several Cartesian degrees of freedom. 

Model 3.1. The Hamiltonian of the total system is defined by the Weyl- 
Wigner symbol 

H {n) +N t (3.1) 

on phase space M" +1 x R n+ . The zeroth components (p,q) of p := 
(p,Pi, ■ ■ ■ , p n ) € M™ +1 and q := (q, qi, . . . , q n ) <G R rl+1 make up the phase- 
space coordinates of the white-noise system, the following components those 
of the heat bath. The deterministic part 

n 

ff(")(p,q) :=p 2 /(2m) + (mu: 2 /2)q 2 + ]T {p 2 j /(2m j ) + {m^/2)^ - qf) 

i=i 

(3.2) 
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of (3.1) couples a harmonic oscillator of mass m > and frequency u > 
bilinearly to a heat bath of n independent harmonic oscillators with masses 
rrij > 0, j = 1, . . . , n, and frequencies < u>\ < u>2 < ■ ■ ■ < u n . The random 
part 



N t (p,d) = N t (q) 



(3.3) 



defines a time- dependent local potential on the phase space of the white-noise 
system only. It is a Gaussian random field with zero mean and covariance 
function 



N t (q)N t ,(q') = S(t-t')C(q-q r ). 

Introducing the diagonal (n + 1) x (n + l)-mass matrix 
M := diag(m, mi, m2, ■ • ■ , m n ) > 

and the (n + 1) x (n + l)-matrix of squared frequencies 



(3.4) 



(3.5) 



n 2 := 



- J2 K?U? -Kiwf 



A 



V 







>o, 



(3.6) 



where «j := yjmj/m, j = 1, . . . , n, one can rewrite (3.2) as a non-negative 
quadratic form 



# (n) (P,q) = ip-M^p+iq-M'/W^q. 



(3.7) 



This model for the heat bath is referred to in the literature as the independent 
oscillator model (Ford et al., 1988a, b) or Caldeira-Leggett model (Weiss, 
1993). 

Following the strategy of the preceding section we define the averaged 
dynamics of the total system as 



T t {n) := lim cxp j- J ds [H( n ) + N ( S T \ .] J , 



(3.8) 



where iV t ^ is a Gaussian coloured noise which converges to the Gaussian 
white noise N t of (3.3) in the limit of vanishing correlation time r N, 0, confer 
Subsection 2.1. 

The validity of the following statement depends only on the quadratic form 
(3.7) of The specific forms of the matrices M and fi 2 are irrelevant. 
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Statement 3.2. The averaged dynamics is a semigroup and admits the 
representation 



T (n) = e _ t£ (») 



cxpj-^ 2 ^ ds (C(0) - C(ift#j n) )) j . (3.9) 



_ffere we /iaue introduced the differential operator 

:= (V 1 / 2 cos^M 1 ^ ^ _ M -i/2 £^f) M -i/2 ^ (3 1Q) 

as i/ie zeroth component of the vector on the right-hand side of (3.10). The 
superscript n accounts for the fact that the Liouvillian which generates the 
reversible dynamics of the deterministic part of the total Hamiltonian is 

C (n) ._ [_ ff (n) ) 

Proof. The semigroup property is established as in the proof of Statement 
2.2. To arrive at the representation (3.9) we remark that the classical phase- 
space trajectories generated by are given by 

e t£<n> p = M 1 / 2 cos(ffi)M- 1 / 2 p - M^ 2 Q sin^M 1 / 2 q , 
e t£(n) q= M-V2 M-^p + M-^cos^M^q 

and that the differential operators IC S and X s introduced in Statement 2.5 
are to be replaced in the multi-dimensional case by vector-valued differential 
operators whose j-th components, j € {0, 1, . . . , n}, read 

{K.)i := (d gj e- sC p) ■ d p + (d q] e- sC <i) ■ 5 q , (3.12a) 
{X B )j := {d Pj c- sC p) ■ d p + {d P] e~ sC q) • 5 q . (3.12b) 

In this notation Definition (3.10) appears as X^ = (X s ) - ■ 

To describe the averaged dynamics of observables we introduce the adjoint 
(T t ( " ) ) t of the semigroup with respect to the standard scalar product 

(f,9):=[ d"+V" + Yr(p,q).9(p,q) (3.13) 

JR"+ 1 xR™+ 1 

for square-integrable functions f,g on phase space R n+1 x M™ +1 . 

Since we are mostly interested in the behaviour of the mean energy, it 
suffices to concentrate on the averaged dynamics of the system observables 
p 2 and q 2 , and thus on the spreading of an averaged state of the white-noise 
system. Again, explicit results on the averaged dynamics of polynomial ob- 
servables may be obtained from a Taylor expansion of the adjoint semigroup. 
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Statement 3.3. 




(3.14a) 




sin(Qs) 



M 



1/2 



) 



00 



} 3 

3.14b) 



3.2 Elimination of the bath variables and the quantum 
Langevin equation 

The second step of the programme outlined at the beginning of this section 
is to eliminate the degrees of freedom of the heat bath assuming that the 
bath is in a thermal equilibrium state. This is achieved with the help of the 
partial expectation value 



of an observable a n of the total system. Consequently, we do not allow for 
initial correlations between the white-noise system and the heat bath. The 
phase-space function in the second line of (3.15) is the Wigner density of the 
canonical equilibrium state of the n bath oscillators at inverse temperature 
P < oo. Note that quantum fluctuations of the heat bath are fully taken into 
account through effectively increased temperatures, defined by the function 



Starting from the Heisenberg equations of motion for the n+ 1 components of 
position and momentum of the total system, one deduces that e* £(n) q satisfies 
a linear quantum Langevin equation (Ford et al., 1965; Benguria and Kac, 
1981; Maassen, 1984; Ford et al., 1987, 1988a; Gardiner, 1991; Jaksic and 
Pillct, 1997a) — derivatives with respect to time are symbolized by superposed 
dots — 





(3.15) 




(3.16) 




(3.17) 



with the friction kernel 



n 




(3.18) 
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and the driving force 

n 

F n (t) := \2{ m j UJ jQj cos(ujjt) + ujjPj sm((jjt)} . (3.19) 
i=i 

The latter behaves with respect to the partial expectation value (•) J a, n as a 
(classical) stationary Gaussian stochastic process with zero mean and cova- 
riance function 



2 



20) 



(F n (t)F n (t')) , n = J2 R ]Js ooB(uj(t f )) =: * Pin (t - f ) . (3. 

Observe that the linear quantum Langevin equation (3.17) in the 
Weyl-Wigner-Moyal representation differs from the corresponding classical 
Langevin equation (Anglin and Habib, 1996; Jaksic and Pillet, 1997b) only 
by the effective temperatures in the covariance function of the driving force. 
Thus we have shown 

Statement 3.4. The classical phase-space trajectories generated by can 
be written as 

e tc(n) p = md t e tc(n) q, (3.21a) 

e tc{n> q=(p/m)G n (t) + qG n (t)+m- 1 / ds G n (t - s) F n (s) . (3.21b) 

Jo 

Here 

f dz e~ izt 
G "W : =L^ <S-z'-izUz),rn > *e*.*>°. (3-22) 

is the retarded Green function (or causal response function) of the quantum 
Langevin equation (3.17) and 

poo . n / ^ \ 

%(z):= dt^ ln {t) = - Y,m^ [—— + ——) (3-23) 
Jo j~i ^ ~ j j ' 

is the Laplace transform of the friction kernel (3.18) defined for z € C , 
Imz > 0. 

Remarks 3.5. i) One deduces for Imz > from the explicit expression 
(3.23) that 

2 2 

Im{ ^ ~ Z } <0<Im{i7»(z)}, (3.24) 

Hence, the denominator of the integrand in (3.22) has no zeros in the 
complex upper half-plane. Moreover, the integrand of (3.22) is analytic 
in this half-plane, and therefore the definition of the Green function G n 
does not depend on e. 



22 



ii) Note that limGJt) = 0, lim G n (t) = 1 and lim G„(i) = 0. 

m ) The Laplace transforms of G n and its derivatives read for z e C , Im z > 
0, 

f°° 1 
G„(z) := / dt c Izt G n (t) = 2 . - , w , (3.25a) 

Jo u 2 - z z -iz~f n (z)/m 

d n (z) = -izG n (z), (3.25b) 
t n {z) = -z 2 G n (z)-l. (3.25c) 

Basically, Statement 3.4 reduces the analysis of the first summands of 
(T^) p 2 and (T t ^)\ 2 , as given by (3.14), to an analysis of the Green 
function G n . The next statement does the same for the second summands in 
(3.14). 

Statement 3.6. Let Ss := {z e C : | Im z\ < 5} be a strip of width 25 > 
centred around the real axis in the complex plane and let f : Ss — ► C be a 
bounded analytic function on Ss with reflection symmetry f(z) = f(—z) for 
all z G Ss- Then, for any < e < S, the 00-matrix element (/(Q)) 00 , where 
f2 stands for the positive square root of the (n+ 1) x (n+ 1) -matrix of squared 
frequencies (3.6), may be expressed in terms of a principal-value integral 

(m) Q0 = - PV / dzzf(z)G n (z). (3.26) 

Proof. We closely follow Ullersma (1966) and Ford et al. (1988b). Using the 
spectral representation Q, = U diag(f2o, • • ■ , ri„)U t , where U is an orthogonal 
(n + 1) x (n + 1)-Matrix and < fio < fl\ < . . . < 0„ are the eigenvalues of 
fl, the left-hand side of (3.26) is rewritten as 

n 

(/(tt)) 00 = E( U °fc) 2 /(^)- ( 3 - 2? ) 

fc=0 

The proof of the following claim is contained in Ford et al. (1988b). 
The function 

e~a(fl:=<Wfl= 2 , .1^, (3-28) 

w -£-iV£7n(v£)/ m 

extends to an analytic function onC\ {Oq, . . . , Moreover, ttj ^ /or 
j 7^ fc, and a has a simple pole at £ = 0^ wif/i residue ResQ2{a} = — (Uofc) 2 
for k = 0, . . . , n. 



Due to the symmetry /(z) = f(—z) there exists an analytic function g on 
S 2 := {£ e C : £ = z 2 ,z e S s } such that g{z 2 ) = f(z). Hence, the 
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(b) 
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Im z 
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K 2 












K 3 




i -n n 


■ -n 


n 




? Re z 







Fig. 2. The poles of the integrand and the integration contour for (a) the integral (3.29) 
and (b) the integral (3.30). 
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above claim and (3.27) lead in combination with the residue theorem to the 
representation 

(/(fi)) oo = i// ffK)a( °- (3 - 29) 

As is depicted in Figure 2(a), the integration contour K encloses clockwise 
all poles of a. It intersects the real axis at £ = — e 2 and at £ = £ 2 , where 
< e < S and £ > Q n . Next we introduce z := y/l; as the square root of £ with 
branch cut along the positive real half-axis and perform a change-of- variables 
in the integral (3.29) 

(/(Q)) =1/ dzzf(z)a(z 2 ), (3.30) 

71"! JK 1 +K 2 +K 3 

see Figure 2(b). In the limit i — > oo the contributions of the vertical line 
segments K\ and K 3 to the integral (3.30) vanish due to the boundedness of 
/, whereas the contribution of the horizontal line segment K 2 converges to 
the principal- value integral in (3.26). ■ 

3.3 The macroscopic limit and the reduced averaged dy- 
namics 

In order to enable the white-noise system to dissipate energy into the heat 
bath, it is necessary to take the macroscopic limit of the heat bath such that 
a continuum of frequencies arises. More precisely, guided by the relation 



7> 



.(*)= [ dv J„Me irt , (3.31) 
Jr 



where J n (y) := | 7] mjUi 2 \8{y — u>j) + 5(y + Wj)], we require that with 

3=1 

the number n of bath oscillators going to infinity, the limit of the friction 
kernel -f(t) := lim 7«(i) exists and admits a representation as the Fourier 

transform of the so-called spectral density J of the macroscopic heat bath 



7(t) - / dv J{u) e wt . (3.32) 

The spectral density describes the strength of the coupling between the sys- 
tem and the heat bath. In what follows, it is supposed to satisfy 

Assumptions 3.7. i) J is continuous, even and strictly positive, J(v) = 
J(-v) > for all v e R . 

ii) There are two constants < V < t/2 and 5 > such that J ad- 
mits an analytic continuation to the subset Tt^, s :— {z e C : Imz < 
and arg(z + i5) ^ [— tt + tp, — tp] } of the complex lower half -plane, see 
Figure 3. 
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Fig. 3. The subset s := {z G C :lmz<0, arg(z + i<5) £ [-tt + V, -V"] } of the complex 
plane. 

m ) J has the following decay properties: 

sup |z 2 J(2;)| < oo and the limit L := lim v 2 J{v), taken along the 

real axis, exists and is approached sufficiently fast, 

lim {\v\ r (v 2 J{v) -L)} = for some r>0. (3.33) 

Note that Assumptions 3.7. i) and 3.7. Hi) imply the integrability of J over 
R. 

Assumptions 3.7 are obviously satisfied for the Gaussian spectral density 
J(v) = Jo expj — ^ 2 /(2o)q)}, Jo, uj n > 0. Another admissible choice for J is 
the Cauchy-Lorentz spectral density J(v) = JoWq/(wq + v 2 ) of the so-called 
Drude model. Both densities approximate Ohmic damping, which is defined 
by J(v) = Jo for all v, the better the larger loq is. Ohmic damping itself 
is excluded by the decay assumption 3.7.m). It leads for ft 7^ to some 
unphysical ultraviolet divergencies, see for example (3.34) and (3.53a) below. 
The reader should be aware that the recent attempt of Ford and O'Conncll 
(1996) to interpret some of these divergencies in a naive distributional sense 
has failed: the proposed distributional derivative of the hyperbolic cotangent 
makes only sense for test functions which vanish at the origin. But this 
implies that the distributional derivative does not differ from the classical 
one. 

For later purpose we state that the covariance function (3.20) of the driving 
force F n in the quantum Langevin equation evaluates in the macroscopic limit 
to 

•Wt) := lim $g „(t) = (ft/2) / d^ J(u) v coth(/3ft^/2) c [l/t . (3.34) 

Accordingly, its Fourier transform differs from that of the friction kernel sim- 
ply by the factor of l//3 e ff(ft^), a fact which reflects the fluctuation-dissipation 
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theorem (Kubo et al., 1991). Furthermore we use the notation 

G(t) := lim G n {t) (3.35) 

n— >oc 

for the macroscopic limit of the Green function of the quantum Langevin 
equation. It does not only depend on the spectral density J of the heat bath, 
but also on the mass m and the frequency u) of the system Hamiltonian. 

After these preparations and recalling the partial expectation value (3.15) 
with respect to the canonical equilibrium state at inverse temperature (3 of 
the heat bath, the reduced averaged dynamics of system observables on phase 
space R x R is defined as 

(^' J) f ■■= ^J( T t n) f •>,,„. ( 3 - 36 ) 

confer (3.8) and (3.15). Its adjoint iz\^' j " 1 with respect to the standard scalar 
product for square-integrable functions on phase space R x R describes the 
reduced averaged dynamics for the states of the white-noise system. 

The starting point for deriving the long-time spreading of states under this 
dynamics is 

Statement 3.8. 

(n ( f' J) y P 2 = ( P G(t) +mqG(t)) 2 + / ds / ds' G(s)$p(s- s')G(s') 

Jo Jo 

+ (-d 2 C)(0) fds (G(s)) 2 , 
Jo 

(3.37a) 

(Tlf' J) ) ] q 2 = (-G(t)+qG(t)Y +m~ 2 fds fds' G{s) $ (s - s') G{s') 
^ m ' Jo Jo 

+ m- 2 (-d 2 C)(0) fds (G(s)) 2 . 
Jo 

(3.37b) 

Proof. The terms in the first line of (3.37a) and (3.37b), respectively, follow 
from taking the partial expectation value (3.15) of the first term in (3.14a) 
and (3.14b), respectively. It is evaluated in the macroscopic limit by using 
Statement 3.4, as well as (3.20), (3.34) and (3.35). The term in the second 
line of (3.37b) emerges in the macroscopic limit from the second term in 
(3.14b) due to 

/sing£) \ = _f dz _ e . =Gn{sh 

V " J 00 JR+ie 2lT 

Here we have applied Statement 3.6 and used the causality condition 

G n (s) = for all s<0, (3.39) 

which is a consequence of Cauchy's integral theorem for analytic functions 
and Definition (3.22). As to the term in the second line of (3.37a), we simply 
refer to (3.38) and to cos(Qs) = d s sin(Qs)/f2. ■ 
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3.4 The long-time limit 

The dissipative character of the reduced averaged dynamics will become ob- 
vious from the long-time behaviour of the Green function G and of the co- 
variance function <Pp. In fact, under our Assumptions 3.7 these functions 
exhibit an exponential relaxation. 

Statement 3.9. There is a decay constant r\ > such that the long-time 
behaviour of the macroscopic Green function (3.35) and of its derivatives is 
given by 

fMo))^^), if u = o t 

\0(c-^), if u>0, 

G(t) = 0(0-"*), (3.41) 

G(t) = ©(e-" 4 ) . (3.42) 

The covariance function <&p from (3.34) is Lebesgue integrable over R . It 
decays exponentially for \t\ — ► oo. For t — ► it either remains finite or 
diverges logarithmically. 

Proof. The Assumptions 3.7 on the spectral density J translate into analyt- 
icity and other properties of j. We summarize them in an 

Assertion. There are two constants < ip < ir/2 and < r\ < 5 such that 7 
can be continued analytically from the upper complex half-plane to the sector 
:= H^ >v U{zeC : Imz > 0}. Moreover, 

sup |27(,z)| < 00, (3.43) 
and the complex function 

z J 2 — z 2 — iz A f(z)/m (3.44) 

has no zeros in H^, v — except for the case u> = in which z — is a simple 
zero. 

Before we show the validity of the Assertion, we demonstrate how it can be 
exploited to prove the Statement. Consider the inverse Laplace transforma- 
tion of the functions G(z) = —izG(z), G(z) = —z 2 G(z) — 1 and G(z) (in the 
case ui > 0), respectively G(z) + (izj(0)/rn) (in the case u = 0). For the 
last function to make sense, observe that 7(0) = ttJ(0) > due to Assump- 
tion 3.7. i). According to the Assertion each of these functions is analytic in 
H^ lV , remains bounded for z — ► — \r} and decays like C(|z _1 |) for \z\ — > 00, 
z e H^, v . Therefore Cauchy's theorem allows to deform the integration 
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contour in the inverse Laplace transformation from the upper complex half- 
plane to the boundary &H^^. But now we are in the situation of the Abelian 
theorem in Doetsch (1950, Chap. 15, §4, Thm. 2), and (3.40)-(3.42) follow. 
To prove the above assertion observe that 

j(z):=if dz'^\, «e^,„ < r? < 5, (3.45) 

coincides in the upper half-plane Im z > with the former definition (see 
(3.23) in the macroscopic limit) 

7(z) = [°°dt c izt 7 (t) = i f dv . (3.46) 

Jo Jr z - v 

This is due to Assumptions 3.7. ii) and Hi) and Cauchy's theorem. Moreover, 
for zo G 1~Lxf),ri the difference quotient (7(2) — 7(^0)) / ( z ~ z o) exists in the limit 
z — > z by Assumption 3.7. m) and dominated convergence. Hence, (3.45) 
provides the desired analytic continuation of 7. The decay property (3.43) 
is evident from (3.45) and Assumption 3.7. Hi). A necessary condition for a 
zero of (3.44) — except for the zero z = in the case u> = — is Im z < as 
follows from Remark 3. 5.2), Re j(y) = ttJ(v) > for eel and Assumption 
3.7.*). Furthermore, (3.43) implies that all zeros lie in a bounded region of 
. Hence, there exist only finitely many of them, because otherwise they 
would have a limit point, and by the identity theorem for analytic functions 
the function (3.44) would have to vanish for all z e H^^- One can therefore 
choose if) and r\ sufficiently small such that there are no zeros at all in H^,. v . 
This completes the proof of the Assertion. 

It remains to prove the properties of the covariance function (pp. Because 

of 

$p(-t) = $/}(t) = h Re |^ Av J W) v coth(/3/uV2) c-' wt ^ (3.47) 

we assume t > from now on. As a consequence of (3.33), J(v) either decays 
like v~ 2 for v — > 00 or it is of the order 0(y~ 2 ~ r ) for some r > 0. In the latter 
case the integral (3.47) exists as a Lebesgue integral for all t > rendering 
<Pp a bounded function. In the former case the integral (3.47) exists as an 
improper Ricmann integral for t > 0, and <Pp has a singularity at t = 0. In 
order to characterize this singularity, fix some s > 0, set L := lim [y 1 J(y)) 

and make the decomposition 

J{v) v coth(/3^/2) =: ^ - ^ c- sv + f(u) (3.48) 

which defines the Lebesgue integrable function / over [0, oof. Hence, / does 
not contribute to the singularity of at t = 0. The singularity of ^ comes 
solely from 

r*,(L-±e-Ae-" t (3.49) 
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and the Abclian theorem in Doetsch (1950, Chap. 13, §3, Thm. 1) reveals 
that it is logarithmic. 

In order to derive the exponential decay of <Pp for t —> oo we proceed as 
in the case of the Green function. Analyticity and decay of J in imply 
together with Cauchy's theorem the representation 

0p(t) = (h/2) [ dv J{v) v coth(/3^/2) e~ il/t 

= (h/2) f dz J(z) z coth(/3fiz/2) e' izt , (3.50) 

which is valid for t > and < rj < min{<5, 2ir/ (f3ti)} ■ The Abelian theorem 
(Doetsch, 1950, Chap. 15, §4, Thm. 2) then yields the exponential decay of 
<£/3 for t — > oo. Thus <Pp is also seen to be Lebesgue integrable over R . ■ 

To formulate the main result of this section, it is illuminating to study first 
the expectation value 

((a))p,j := lim («;<">, a) (3.51) 

n — >oo 

of a system observable a with respect to the thermal equilibrium state of the 
deterministic part of the total Hamiltonian in the macroscopic limit of 
the heat bath. The scalar product in (3.51) was defined in (3.13). Recalling 
the effective temperatures (3.16), the Wigner density of the corresponding 
pre-limit state may be written as 

»<">(p,q) - dC %^ll m)) exp{-| P • M-V^ eff (^)M-V 2p 

- \ q • M 1 / 2 n 2 /3 eff (^0)M 1 / 2 q} . (3.52) 

Statement 3.10 below compiles three equivalent expressions for the thermal 
expectation values ((p 1 ))/3,j and (((j' 2 ))/3 i j, which are available in the literature, 
see for example Haake and Reibold (1985), Ford ct al. (1988b), Grabert et 
al. (1988, Table 2) and Weiss (1993, Chap. 6). Nevertheless we will prove 
them, since Statement 3.6 provides a quick way to do so. 

Statement 3.10. The following identities hold under the Assumptions 3.7 
on the spectral density: 

iP 2 h,J = % I & ■/(") ^ \G{vf coth(/?fo//2) 

= ^? pv[ dv (is 2 G(v) + 1) cath(/?fa//2) 

+ 00 

= ? E (i-^ahi)) 

f l=-oo 



(3.53a) 
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and, ifu>>0, 

({q 2 h,j = ^ J/" J(v)v\G(v)\ 2 coth(/?fo//2) 
pv[ Av G{v) coth(/3/M//2) 



ft 



27rim 

+00 

= ^3 E ^hl)- (3.53b) 

i— — 00 

ffere, z/; := 2irl/((3h) are the even Matsubara frequencies, G(v) denotes the 
analytic continuation of G(z) to the real axis, and the principal-value pre- 
scription in the middle lines of (3.53a) and (3.53b) has to be applied at the 
origin. 

Remark 3.11. As is indicated by the notation, the expectation values ((p 2 ))/3,j 
and ((<7 2 ))/3,,/ depend on the inverse temperature (3 and on the spectral density 
J of the heat bath. Moreover, they do also depend on Planck's constant H. 
From the series representations in (3.53) one can derive the following chain 
of inequalities 

r 1 (uj/ujj) 2 } mui 2 „ ,.. 1 



2/3' 2f3 eS (hwj) i ~ 2 xxy " P ' J ~ 2/? eff (M 



^L^^wk^r (3 ' 54) 



where wj := yj lo 2 + 7(0)/m and 7(0) = f R dv J{v). As a consequence, one 
can easily deduce the values of ((p 2 }} / s ! j/2m and raw 2 (( q 2 ))p t j/2 in the low- 
friction limit J \ or the classical limit h \ 0. For J \ both expectation val- 
ues tend to the quantum-mechanical canonical equilibrium value l/2[3 c s(ftiu) 
with respect to the harmonic-oscillator Hamiltonian p 2 /(2m) + (moo 2 /2)q 2 . 
For H\0 they tend to 1/2/3 for all admissible J, as dictated by the equiparti- 
tion theorem. The latter expression emerges also in the quantum-mechanical 
case h 7^ as an asymptotic value for high temperatures f3\0. 

Proof of Statement 3.10. In order to calculate ((p 2 ))/3,j one has to consider 
the macroscopic limit of the expectation value 



(„,("), p 2 ) = ^ (ttcoth(/?MV2)) 00 . (3.55) 



Before doing this we rewrite the right-hand side of (3.55) with the help of 
Statement 3.6. To this end, observe that the function 

z ~ /(*) = i + x-2 z 2 coth(/?ftz/2) (3.56) 
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fulfills the assumptions of Statement 3.6 for all < e < 2ir/ ((3ft) < A. Hence, 

(«,(») ,p2) = ^ lim f dz — z 2 G n (z) coth((3hz/2) 
2tti a^oo J R+i£ 1 + A l z z 

= S / dZ (^"^ + X ) coth (/3^/2) 

-^lim / dz ^ , (3.57) 

27Ti A^oo J R+i£ 1 + A- 2 Z 2 V ' 

where we have taken advantage of the fact that z 2 G n (z) + 1 falls of like |z|~ 2 



f , coth(/3ftz/2) /■ , coth(/37ii//2) . ^ r , ,„ 

/ d2 u - 2 2 = pv / d ^ i T A -2,> - 7r 1 Rcs z=0 {coth(/3ftz/2)} 

•/R+ie 1 + A Z J R 1 + A f 

= -2iri/{/3ft) (3.58) 

we conclude 

((P 2 ))/3,J = tj-t / dz (z 2 G(z) + 1) coth(/3fe/2) + 27ri/C9ft) 

Z7T1 iJR+ie 

= !r? pv / d ^ (^ 2 G(^) + 1) cath(/?fe//2) . (3.59) 

The principal-value prescription has to be taken at the origin, and G(v) 
denotes the analytic continuation of G(z) to the real axis. This continua- 
tion exists due to Assumptions 3.7. The integral representation in the top 
line of (3.53a) follows from (3.59) by using the antisymmetry of the hy- 
perbolic cotangent, the reality condition G(— v) = (G(v)) and lmG(u) = 
7T77J -1 J(v)v\G(v)\ 2 . The series representation in the bottom line of (3.53a) 
follows from evaluating (3.59) with the help of the residue theorem and the 
Mittag-LefHer expansion of the hyperbolic cotangent 

coth(/3fe//2) = - ]T (3.60) 

l— — oo ^ 

Here, vi := 2irl/(f3H) are the even Matsubara frequencies. Note that the pole 
for I — 0, which lies on the integration contour, counts only half due to the 
principal-value prescription. 

In order for ((g 2 ))^j to exist we assume to > and start out with 

We pick three constants 0<e<A<e< 2ir/(pft) and apply Statement 3.6 
to the function 

z " f( z ) = ~vT~ — 2 coth(/3fiz/2) (3.62) 
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which gives 

(w (n \q 2 ) = JL- Inn J 6z ^ G n{ z) coth^/2) 



L+ie 



lim <^ / dz J G n {z) coth(/3fe/2) 



27rim axo [7 R+ii ? A 2 + z 2 

+ 7riAG„(iA) cot(/3fiA/2) 



27rim 



( / dz G„(z) coth(/3ftz/2) + ^ G„(0)1 . (3.63) 



The second equality in (3.63) follows from deforming the integration contour 
across the pole at z = iA. Performing the macroscopic limit n —* oo and 
subsequently the limit e\ 0, one arrives at 

((9 2 »/3,/ = 5-?- pv/ div G(i/) coth(/?fo//2) . (3.64) 

The integral representation in the top line of (3.53b) and the series represen- 
tation in the bottom line of (3.53b) are derived from (3.64) in the same way 
as was demonstrated for {(p 2 }}f3,j- ■ 

As the main result of this section, we describe the spreading of states under 
the reduced averaged dynamics of Model 3.1 in the long-time limit. 

Statement 3.12. If the spectral density J of the macroscopic heat bath sat- 
isfies Assumptions 3.7, then the second moment of the particle momentum 
approaches a finite value under the reduced averaged dynamics (3.36) of Model 
3.1 

lim (tt^jV = ((p'h.j + (-#7)(0) [ P v 2 \G{ V )\ 2 • (3.65) 

In the harmonically bound case u > 0, the same is true for the second moment 
of position 



lim (nr ] )\ 2 = U%, + / P (3-66) 



t^oo v 1 ' " IJ, ° ' to 2 J K 2tt 

In the unbound case w = the limit (3.66) is to be replaced by the diffusive 
behaviour 



& ^~ L ~t~ = urn v 1 + ) ■ (3 - 67) 

Above we have used the thermal expectation values (3.53). The results imply 
in particular that the mean energy {TZ ( f' J ^{p 2 /{2m) + (mw 2 /2)g 2 ) , respec- 
tively (lZ^f ,J ^ p 2 / (2m) in the unbound case uj — 0, remains finite in the 
long-time limit. 
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Remarks 3.13. i) The right-hand sides of (3.65) - (3.67) are constant func- 
tions on phase space 1x1. Consequently, the asymptotic spreading 
of states is universal in the sense that it does not depend on the initial 
state. 

ii) The terms which are due to the white noise in (3.65) - (3.67) depend 
only on the curvature of its covariance function at the origin and, via 
the Green function G, on the spectral density J of the heat bath. On 
the other hand, the thermal equilibrium expectation values (( p 2 ))/3,j 
and ((q 2 }}/3,j inhibit also a dependence on Planck's constant. But it is 
only in the subsequent low-friction limit J \ that they reduce to the 
quantum-mechanical canonical equilibrium values with respect to the 
system Hamiltonian alone, see Remark 3.11 and confer also Ford et al. 
(1965), Davies (1973), Benguria and Kac (1981) and Maassen (1984). 

Hi) In Jayannavar and Kumar (1982) and Heinrichs (1983a) it has already 
been conjectured that the interplay of white noise and dissipation would 
lead to a diffusive spreading in position if u = 0. Eq. (3.67) provides 
the corresponding quantitative statement. It is to be compared with 
the t 3 -behaviour (2.60e) of T^q 2 in the situation without dissipation. 
The ^.-independence of the diffusion constant on the right-hand side of 
(3.67) reveals the classical character of this phenomenon. 

iv) In the absence of white noise, C(0) = 0, Eqs. (3.65) - (3.67) reduce 
to known results for the damped harmonic oscillator (lu > 0) and the 
damped free particle (ui = 0), respectively. In particluar, the diffusion 
constant on the right-hand side of (3.67) is then given by an Einstein 
type of relation (Sutherland, 1905; Einstein, 1905; Kubo et al., 1991). 
In the above statement we have given precise and explicit assumptions 
on the spectral density under which these results are valid. At zero 
temperature, (3 = oo, or for a spectral density with J(0) = or J(0) = 
oo, which is excluded by our assumptions, non-exponential relaxation 
and anomalous diffusion is expected, see for example Haake and Rcibold 
(1985), Grabert et al. (1988) and Weiss (1993, Chap. 6). 

Proof of Statement 3.12. Statement 3.9 implies that the first term on the 
right-hand side of (3.37a) and (3.37b), respectively, does not contribute to 
the leading behaviour in the limit t —* oo. For the case u) = the leading 
behaviour of (3.37b) is determined by the limit (3.40). The result (3.67) 
follows from it with the help of the relation 




(3.68) 



The integral J *ds (G(s)) 2 in (3.37a) converges for t — > oo. For w > 
the same is true for the integral f*ds (G(s)) 2 in (3.37b). Using Parseval's 
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formula we thus conclude 
dv 

./ o 27T 



ds (GM) 2 = / ^^ 2 |GM| 2 (3.69) 



and, if u) > 0, 

j™ds(G(s)) 2 = j^\G{v)\\ (3.70) 

This explains the last term in (3.65) and (3.66), respectively. Finally, the 
convolution theorem, Definition (3.34) and Statement 3.10 imply 

r-OO r-OC 

ds &s'G(s)$f}(s- s')G(s') 
Jo Jo 

= ? / di> J(v) \G{v)\ 2 coth(/3^/2) 

= (P 2 h,i (3-71) 

and, if uj > 0, 

fOO /'OO 



™- 2 



/•OO /• OO 

/ ds / d8'G(8)$fi(8-sf)G(8') 

Jo Jo 

[ dv J(v) v |GV)| 2 coth(/9fii//2) 



_ 2m" 2 

= U 2 h,J- (3-72) 



4 Outlook 

The present paper is devoted to the averaged dynamics of white-noise sys- 
tems. Natural extensions of this work concern the fluctuations around the 
averaged dynamics and the effects of a non-vanishing correlation time, that 
is, coloured-noise perturbations. Both problems have not been treated satis- 
factorily so far. Diverse attempts, both analytical (but based on uncontrolled 
approximations) and numerical, can be found in the literature. To our knowl- 
edge, all of them treat only the case without dissipation. Unfortunately, their 
results are not free from contradictions. It may therefore be of help to compile 
some of these references and to relate their findings. 

In the preceding sections we analyzed the averaged dynamics w(p, q) i— ► 
wt (p, q) in order to study the average (wt , a) of the quantum- mechanical 
expectation value (w*,a) of an observable a at time t. Fluctuations around 

this average are described by the higher moments (w t7 a) N , N > 2, for which 
Jensen's inequality 

(w t ,a) N >(jM) N (4.1) 
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provides a (generally rather poor) lower bound. The appropriate tool for 

studying the iV-th moment (wt,a) N is the averaged dynamics of a non- 
interacting iV-particle white-noise system. In particular, it provides the 
relevant information about the mapping w(pi,qi) ■ ... ■ w(pn,<1n) * 
Wt(pi,qi) ■ .. . ■ Wt(pN ,Qn)- This dynamics is again given by a semigroup, 
whose generator can be constructed explicitly. However, since one has not 
succeeded in representing the semigroup as in Statement 2.5, which was the 
starting point for all the exact calculations, one had to resort to approx- 
imations. Stintzing (1992) has carried out an analysis up to second-order 
perturbation theory for N — 2 particles on a lattice with nearest-neighbour 
hopping and a site-diagonal white-noise perturbation. Taking the lattice po- 
sition q — (qi, . . . , qd) as the observable, he finds an asymptotic result which, 
in our notation, would read 

(w t ,q) 2 ~t° for t^oo. (4.2) 

This result is shown to be independent of the number d of space dimensions 
of the lattice. By appealing to an analogy to a classical master equation, he 
even claims that (4.2) represents the exact asymptotic behaviour. However, 
(4.2) should be contrasted with the numerical calculations of Bouchaud et al. 
(1992) and Saul et al. (1992) which suggest a result different from (4.2) 

M)^{l t ; d d ZY for t^oo. (4.3) 

Other aspects of the averaged dynamics of many-particle white-noise systems 
can be found in Girvin and Mahan (1979) and Jauho (1987). 

Coloured-noise perturbations have attracted much attention in the last 
decade. On the one hand they give rise to less idealized models, on the other 
hand they provide a variety of interesting features due to the existence of 
an additional time scale, the correlation time r of the noise. Because of 
this complexity, there are only few exact analytical results available. More- 
over, these are limited to models in which the noise has some additional 
nice properties. We only mention a Markovian noise in time (Cheremshant- 
sev, 1982; Pillet, 1985, 1986; Tcheremchantsev, 1997), rank-one perturbations 
with coloured- noise coupling coefficients (Chvosta, 1991; Neu and Speicher, 
1994) and coloured-noise couplings to some deterministic external force (Ki- 
tahara and Haus, 1979; Inaba, 1981). For classical systems, however, there 
exists a more detailed picture, see for example the review by Hanggi and 
Jung (1995). Quantum-mechanical models similar to ours with a Gaus- 
sian coloured-noise have been studied by Golubovic et al. (1991), Rosenbluth 
(1992), Hcinrichs (1992), Lebedev et al. (1995) and Hcinrichs (1996). The 
starting point for the calculations in Golubovic et al. (1991) is the corre- 
sponding classical model which is analyzed within the Martin-Siggia-Rose 
formalism. After performing a short-time expansion they derive the asymp- 
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totic long-time behaviour 



{w t 



V 2 ) ~ { 



(wt,q 2 ) 



t 2 ^ , d=l 
t 1 ' 2 , d>2 

for t — > oo , (4.4) 

i 12 / 5 , d=l 
t 9 / 4 , d>2 



which, according to a quasi-classical argument, should not be modified by 
quantum fluctuations. Golubovic et al. (1991) also perform a numerical sim- 
ulation in order to support (4.4). In a Comment, Rosenbluth (1992) criticizes 
the derivation of (4.4), as it does not guarantee the conservation of proba- 
bility. His alternative derivation, which is not presented in detail, leads to a 
modification of (4.4) in more than one space dimension 



(w t ,p 2 ) ~ t 2 ' b 

for t — > oo . (4.5) 

7 =n ft 12/5 , d=l 

A genuine quantum-mechanical treatment of the model is given by Hcinrichs 
(1992). He succeeds in averaging the von Neumann equation up to first order 
in r/t yielding 



(w t ,p 2 )~t {l + 0((r/t) 2 )} 



for t -» oo . (4.6) 



(^.^-^{l + f (r/t) + 0((r/t) 2 )} 



Using a similar method, Heinrichs (1996) also investigates the corresponding 
behaviour at intermediate times. Lebedev et al. (1995) argue on rather spec- 
ulative grounds that (4.4) should be correct up to very long times, but in the 
limit t — > oo the exponents should be replaced by smaller ones. 

We conclude that the effects of a non-vanishing correlation time are con- 
troversially discussed at present. In order to arrive at a unifying picture, 
it seems that suitable bounds and other controlled approximations to the 
averaged dynamics are needed. 



Appendix A The Weyl-Wigner-Moyal representation 

In this Appendix we intend to give a brief overview of the linear phase-space 
representation of quantum mechanics which dates back to Weyl (1928), Wig- 
ner (1932) and Moyal (1949). More detailed information and proofs of the 
results mentioned here can be found in the review articles Tatarskii (1983), 
Balasz and Jennings (1984) and Hillary et al. (1984). Extensions of the for- 
malism in order to describe many-particle systems are dealt with in Gadella 
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and Nieto (1993, 1994), extensions to a non- Euclidean phase space in Bayen 
et al. (1978) and Kasperkovitz and Peev (1994). Some more mathematical 
issues are contained in Pool (1966), Daubechies (1980, 1983), Folland (1989), 
Lieb (1990) and Fedosov (1996). The articles Molzahn and Osborn (1994), 
Osborn and Molzahn (1995), Arai (1995) and Robinson (1993) are concerned 
with quasi-classical properties and the classical limit. 

For simplicity we will restrict ourselves to describe a single spinless quan- 
tum particle which has the Euclidean line R as its configuration space. The 
generalization to several Cartesian degrees of freedom is merely a matter of 
notation. In the standard Hilbcrt-space formulation states and observables of 
the particle are represented by certain linear operators on L 2 (R ) , the Hilbert 
space of Lebesgue square- integrable complex- valued functions on R . The link 
to phase space is provided by 

Definition A.l. Let F be a linear operator on L 2 (R). The mapping W : 
F h-> W(F) := /, 

/:lxR-»C, f(p, q) := f dr c [pr / h (q - r/2\F\q + r/2) (Al) 

is called Weyl-Wigner mapping. The phase-space function f is called the 
Weyl-Wigner symbol of the operator F. In (Al) we have used the Dirac 
notation for the integral kernel of an operator in the position representation. 

Pool (1966) shows that the Weyl-Wigner mapping induces an isometric iso- 
morphism between the space of Hilbert-Schmidt operators on L 2 (R ) and the 
Hilbert space L 2 (R x R ) of square-integrable complex- valued functions on 
phase space. For an extension of (Al) to the Banach space of bounded oper- 
ators on L 2 (R), see Daubechies (1980). For the purpose of the present paper 
we may safely ignore these more subtle issues, and we will use (Al) for all 
operators F for which it makes sense, possibly in a distributional one. 

Some important properties of the Weyl-Wigner mapping are listed without 
proof in 

Statement A. 2. i) Compatibility with the complex-linear structure: 

W(aF + [SG^) = aW(F) + /3W(G)* = af + (3g* , (A2) 

where a, [3 £ C, is the Hilbert adjoint of G on L 2 (R) and g* is the 
complex conjugate of g. 

ii ) Operator products correspond to star products of symbols ( the arrows 
indicate on which symbol a particular differentiation is meant to act): 

W(FG) = /exp|-| (d p d q -d q d p )\ g=:f*g. (A3) 
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Hi) Isometry: 

Tracc(FtG) = (2nh)- 1 f dpdq f*(p,q) g(p,q) =: (2^fi)- 1 (.f, ff ) . 

JRxR 

(A4) 

iv) Weyl quantization as the inverse mapping: 

F = W~\f)= ( ^l c ^-^/\f(p,Q + q/2) 

= f{-ihd q >,-ihdp>) e i(9 ' p+p ' Q)/ft . (A5) 

p'=0=q' 

Here P and Q denote the momentum and position operator, respec- 
tively. 

v) Diagonal matrix elements by integration: 

i 2% q) = {q ™ ' i 2% q) = {PIFIP) • (A6) 

Remark A. 3. According to Statement A.2.n) the symbol of the standardized 
commutator is given by 

W((i/h)(FG-GF))={i/h){f*g-g*f) 

=:[f,9}. (A7) 



It is called Moyal bracket. If either / or g is a quadratic polynomial in p and 
q, then the Moyal bracket reduces to the Poisson bracket 

[f,g} = (d p f)(d q g)-(d q f)(d p g). (A8) 

Regardless of this assumption on / and g, the equality (A8) is always true in 
the classical limit H\ 0. 

There are uncountably many ways of associating a phase-space function 
with a quantum-mechanical operator and vice versa (Cohen, 1966; Agarwal 
and Wolf, 1970; Fink and Leschke, 1993). The Weyl-Wigncr mapping, re- 
spectively the Weyl quantization, provides the one that treats position and 
momentum in a totally symmetric way. This is illustrated by the following 

Examples A.+. ») VW 1 ^)) - f(Q) , W^/fr)) = f(P), 
ll) W-^pfiq)) = \ {Pf(Q) + f(Q)P), 

ill) W(Pf(Q))=pf(q) - (ifc/2) (d q .f)(q), 



iv 



) W{Pf{Q)P) = p 2 f(q) + (h 2 /4) (d 2 q f)(q). 
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Phase-space functions which represent quantum-mechanical states are called 
Wigner densities. By convention we say that a phase-space function to is a 
Wigncr density if it can be written as 

w = (2Trh)- 1 W(W), (A9) 

where W is some state operator, that is, W = W> > and Trace W = 1. 

Statement A. 5. Wigner densities are bounded, \w(p, q)\ < (irh)^ 1 , normal- 
ized, (w, 1) = 1, and square-integrable, (w,w) < (2wfi)~ 1 . Moreover, one has 
(w,w) = (2-Kh)~ 1 if and only if w represents a pure state. 

In contrast to probability densities, however, Wigner densities may also take 
on negative values. This fact is a consequence of Statements A.2.i) and v), 
see e.g. Wigner (1971). A simple example is provided by the pure state 
W — IV'KV'I whose position amplitude (q\tp) is the positive square root of 
a uniform probability density over some interval. This example shows in 
addition that Wigner densities need not be Lebesgue integrable (Daubechies, 
1983). 

Below we give two intrinsic characterizations of Wigner densities. The first 
one relies on the factorization W — (W 1 ^ 2 ) (Shirokov, 1979), the second 
one reflects the corresponding properties of state operators (Narcowich and 
O'Connell, 1986). 

Statement A. 6. A phase-space function w is a Wigner density, if and only 
if there exists a real-valued phase-space function u with (u, u) — 1 such that 
w = u*u. 

Statement A. 7. A phase-space function w is a Wigner density, if and only 
if it has the following three properties: 

1. w — w* , 

2. (w,l) = l, 

3. (w,v) > for all Wigner densities v. 

Statement A. 6 supports that Wigner densities converge to probability densi- 
ties in the classical limit H \ (at least in some weak sense). From a rigorous 
mathematical point of view, however, such questions seem to be highly non- 
trivial. This is due to severe technical difficulties which arise in part from 
the missing Lebesgue integrability of Wigner densities. Related aspects of 
the classical limit of quantum-mechanical phase-space representations may 
be found in Robinson (1993), Arai (1995) and Werner (1995). 
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